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We present a signature-based search for anomalous production of events containing a photon, two 
jets, of which at least one is identified as originating from a b quark, and missing transverse energy 
(J? T ). The search uses data corresponding to 2.0 fb _1 of integrated luminosity from pp collisions 
at a center-of-mass energy of y/s = 1.96 TeV, collected with the CDF II detector at the Fermilab 
Tevatron. From 6,697,466 events with a photon candidate with transverse energy Et > 25 GeV, we 
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find 617 events with J5 T > 25 GeV and two or more jets with Et > 15 GeV, at least one identified 
as originating from a b quark, versus an expectation of 607 ± 113 events. Increasing the requirement 
on E T to 50 GeV, we find 28 events versus an expectation of 30 ± 11 events. We find no indications 
of non-standard-model phenomena. 

PACS numbers: 13.85Qk, 12.60Jv, 14.80Ly 



I. INTRODUCTION 



Within the standard model of elementary particle physics (SM) there are six flavors of quarks, six flavors of leptons, 
and four vector gauge bosons, with a hierarchy of couplings and masses. The Fermilab Tevatron, with a center-of-mass 
energy of 1.96 TeV, can produce all of the known quarks and vector bosons. Over the course of years of data-taking 
using the CDF detector we have developed a suite of largely data-driven methods by which we estimate the 
efficiencies and backgrounds associated with the identification of charged leptons, heavy flavor quarks (b or c quark), 
electroweak gauge bosons (photon, W ± , and Z°), and the presence of neutrinos, identified generically by missing 
transverse energy (E T ) Q. The ability to identify these 'objects' in events and to estimate their efficiencies and 
backgrounds has led to the development of signature-based searches at the Tevatron, in which one defines a priori the 
objects an event is required to contain, and then compares observations to expectations @, 0, 0, IS IS ■ The model 
tested in these searches is the SM, which is predictive and falsifiable; any deviation from the SM predictions would be 
a signal of new phenomena. The advantage of this strategy is that only once such a signal has been established would 
the investment be made in generating detailed predictions of the many possible models for the new phenomena. 

We describe here a search for new physics in the inclusive 767 channel using 2.0±0.1 fb _1 of integrated luminosity 
at a center-of-mass energy of y/s = 1.96 TeV, collected between February 2002 and May 2007. A similar search was 
originally performed in Run I using « 85 pb _1 of integrated luminosity flOj . Our search in Run II is part of a broad 
effort at CDF to study rare event signatures involving photons for any non-SM sources [E0I1]- The SM processes, 
cither with a radiated photon or where the charged lepton is misidentified as a photon, are expected to contribute 
ss 2% (tt — > Ivjjbb) and < 1% (Wbb — > Ivbb and Zbb — > vvbb) to the measured rate ll|. Because the SM contributions 
to the jbj^T final state are highly suppressed, for an ideal detector the signature provides an excellent place to look 
for new phenomena. In reality, we expect additional events from processes such as 7 + jets and bb production in 
which mismcasuremcnts of the jet energy induce E T . 

The outline of the paper is as follows. Section HT1 briefly describes the CDF II detector. The selection of events with 
photons, jets, jets from a heavy flavor quark (b or c quark), and missing transverse energy is described in Sec. IIIII The 
estimation of backgrounds to the search sample is presented in Sec. IIVI Section [V] describes the sources and estimates 
of systematic uncertainties on the numbers of events from backgrounds. The results of the search, including the effect 
of additional selection criteria, are presented in Sec. IVIl Section IVlII presents the conclusions. 



II. THE CDF II DETECTOR 



The CDF II detector is a cylindrically-symmetric spectrometer designed to study pp collisions at the Fermilab 
Tevatron. The detector has been extensively described in detail in the literature [J. Here we briefly describe the 
detector subsystems relevant for the analysis. 

Tracking systems are used to measure the momenta of charged particles, to reconstruct primary and secondary 
vertices, and to trigger on and identify leptons with large transverse momentum Q. Silicon strip detectors [l2| and 
the central outer tracker (COT) [l3| are contained in a superconducting solenoid that generates a magnetic field of 
1.4 T. The silicon strip system provides up to 8 measurements in the r — <f> and r — z views [2| and covers the track 
reconstruction in the region \r\\ < 2. The COT is an open-cell drift chamber that makes up to 96 measurements along 
the track of each charged particle in the region \r}\ < 1. Sense wires are arranged in 8 alternating axial and ±2° stereo 
super-layers. The resolution in pt, cr PT /pT, is ~ 0.0015 pr ■ GeV -1 • c for tracks with only COT measurements, and 
« 0.0007 • GeV" 1 • c for tracks with both the silicon and COT measurements. 

Calorimeters are segmented with towers arranged in a projective geometry. Each tower consists of an electromagnetic 
and a hadronic compartment [3, HE Gil s covering the central region, \rj\ < 1.1 and the 'end plug' region, 1.1 < |?7| < 
3.6. The central electrogmagnetic calorimeter (CEM) and central hadronic calorimeter (CHA) arc in the central 
region while the plug electromagnetic calorimeter (PEM) and plug hadronic calorimeter (PHA) are in the 'end plug' 
region. In this analysis, a high-E^ photon is required to be identified in the central region, where the CEM has a 
segmentation of 15° in and w 0.1 in 77 [l|, and an E T resolution of a(E T )/E T w 13.5%/V^r/GeV © 2% We 
further require a high-E^ jet to be identified in the central region, where the jet energy resolution is approximately 
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ctrjO.1- Et (GcV) © 1.0 GeV [l7|. Two additional systems in the central region with finer spatial resolution are used 
for photon identification in this analysis. The central strip system, CES, uses a multi-wire proportional chamber with 
1.67- and 2.01-cm-wide cathode strips and a wire spacing of 1.45 cm to make profile measurements of electromagnetic 
showers at a depth of 6 radiation lengths (approximately shower maximum). The central preshower detector, CPR, 
located just outside the solenoid coil on the front face of the CEM, separates single photons from the photon pairs 
from 7T° and rp decays on a statistical basis, as described in Sec. IIV Al In 2005 the CPR was upgraded from the Run 
I configuration of wire proportional chambers, similar to those used in the CES, to a fast scintillator system with a 
segmentation of 12.5 cm in <fi and 12.5 cm in z [l6[. The finer segmentation in z reduces the probability of a random 
hit from the underlying event and multiple interactions by a factor of four, thereby improving the performance of the 
preshower detector in higher luminosity beam conditions. 

Muons are identified using the central muon systems [l8[: CMU and CMP for the pseudorapidity region of \r\\ < 0.6, 
and CMX for the pseudorapidity region of 0.6 < \r)\ < 1.0. The CMU system uses four layers of planar drift chambers 
to detect muons with px > 1.4 GeV/c. The CMP system consists of an additional four layers of planar drift chambers 
located behind 0.6 m of steel outside the magnetic return yoke, and detects muons with pr > 2.2 GeV/c. The CMX 
system detects muons with pr > 1.4 GeV/c with four to eight layers of drift chambers, depending on the direction of 
the muon. 

The luminosity is measured using two sets of gas Cerenkov counters [l9l ]. located in the region 3.7 < \r/\ < 4.7. The 
total uncertainty on the luminosity is estimated to be 5.9%, where 4.4% comes from the acceptance and operation of 
the luminosity monitor and 4.0% from the calculation of the inelastic pp cross-section [20l |. 

A three-level trigger system [2l[ selects events to be recorded for further analysis. The first two trigger levels consist 
of dedicated fast digital electronics analyzing a subset of the full detector data. The third level, applied to the full 
data of those events passing the first two levels, consists of a farm of computers that reconstruct the data and apply 
selection criteria consistent with the subsequent offline event processing. 

III. EVENT SELECTION 

An initial sample of events enhanced with high energy photons is collected using a trigger that requires a high 
energy isolated cluster in the electromagnetic calorimeter [22| . We require events to have a primary vertex with 
\z\ < 60 cm. The offline selection criteria require a central (\rj\ < 1.1) photon with Et > 25 GeV, two jets with 
\rj\ < 2.0 and Et > 15 GcV, at least one of which is identified as originating from a b quark (6-tagged), and missing 
transverse energy greater than 25 GeV, as described in more detail below. The selection is inclusive; i.e. we allow 
extra objects (jets, photons, leptons) in the events. 

The photon is required to satisfy the same identification requirements as in previous CDF high-ET photon analy- 
ses [23| . Namely, the photon candidate is required to have no associated track with pt > 1 GeV, at most one track 
with pt < 1 GeV pointing at the calorimeter cluster [24| , good profiles of electromagnetic e nerg y measured in both 
transverse dimensions at shower maximum, and minimal leakage into the hadron calorimeter [25[ . Photon candidates 
identified via these cuts are referred to as "standard" photons. 

Jets are reconstructed using the jetclu cone algorithm (2|| with cone radius R = yj A0 2 + Arf = 0.4. Starting 
from seed locations corresponding to calorimeter towers with Et > 1 GcV, all nearby towers with Et > 0.1 GeV 
are used to search for stable cones. To resolve ambiguities with overlapping cones, cones sharing an energy fraction 
greater than 0.75 are merged into a single jet; otherwise the shared towers are assigned to the closest jet. We apply 
a jet energy scale (JES) correction [27[ such that the measured Et is on average equal to the summed Et of the 
particles from the pp interaction within the jet cone. At least one of the jets must be 6-tagged using the SECVTX 
algorithm [28[, which searches for displaced vertices using the reconstructed tracks inside the jet cone. 

Missing transverse energy E T is calculated Q from the calorimeter tower energies in the region \r)\ < 3.6. Corrections 
are then applied to the E T for (i) the calorimeter response for identified jets [27| . and (ii) the presence of muons with 
Pt > 20 GeV. We require the corrected E T to be greater than 25 GeV and minimize the number of events with 
mismeasured E T by requiring the difference in azimuthal angle between any jet and the E T , A</>(jet, E T ), to be 
greater than 0.3. The requirement AR > 0.4 is imposed on all combinations of the photon and the two selected jets, 
namely 7J1, 7J2, and jijV One of the two jets is the leading 6-tagged jet, and the other is the next-to-leading 6-tagged 
jet if one exists, or the leading non-6-tagged jet if not. 

Table [J summarizes the event selection. The final 7&?|?t sample with E T > 25 GeV corresponds to one part in 10 4 
of inclusive high- Et photon events. We will refer to this sample of 617 events as the "search" sample. 
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TABLE I: Summary of the event selection. The selection of a central photon includes the requirement of the inclusive photon 
trigger, the selections on the z-vertex, and Et{^) as described in Sec. IIIII The selection A_R > 0.4 is required for each pair of 
TJii 7J2, and j 1 j 2 . 



Cut Events 

Photon with E T > 25 GeV, \n\ < 1.1 6,697,466 

2 Jets with E T > 15 GeV, \r)\ < 2.0 1,944,962 

AR > 0.4 for 7ji, jj 3 , and jij 2 1,941,343 

E T > 25 GeV 35,463 

A(j>(jet, E T ) > 0.3 18,128 

> 1 secvtx 6-tag 617 



IV. BACKGROUND PREDICTIONS 



To understand the composition of the search sample of 617 events, we could in principle use Monte Carlo (MC) sim- 
ulations to estimate the absolute numbers of events with real or misidentified photons, and real or misidentified heavy 
flavor. However, this method would result in large systematic uncertainties on the number of events due to theoretical 
uncertainties on the production cross-sections and difficulties in modeling misidcntifications. Wc have consequently 
developed a data-driven strategy that uses the Monte Carlo judiciously to minimize systematic uncertainties. 

As mentioned in Sec. HI SM processes with final state neutrinos are not expected to contribute significantly to our 
search region. We check this expectation by vetoing events that have any high pt isolated tracks, effectively removing 
any contribution from processes involving leptonic decays of vector bosons. Isolated tracks are defined as tracks with 
Pt > 20 GeV having an isolation fraction larger than 0.9, where the isolation fraction is defined as 

track 

f _ Pt 

J is 



Further details of the isolation calculation are given in Ref. (2j|. After the application of the isolated track veto, the 
observed number of events decreases from 617 to 600. This decrease is consistent with the r* 3% expectation obtained 



from Monte Carlo simulation [11 1 



We define four categories of dominant background events, all of which have missing transverse energy primarily 
arising from energy mismcasurcmcnt. Wc obtain the kinematic shapes and normalizations of each category separately. 
The four categories are a) misidentified photons, referred to as "misidentified 7" , b) true photon plus light quark jet 
misidentified as heavy flavor, referred to as "true 7, misidentified b" , c) true photon plus true b-quark jet, referred to 
as "76" , and d) true photon plus true c-quark jet, referred to as "7c" . 

The misidentified 7 background is estimated from the data sample itself by using cluster-shape variables from the 
CES and hit rates in the CPR. This technique (the CES/CPR method) allows the determination of the number 
of photon candidates in the sample that are actually misidentified jets as well as the corresponding shapes of the 
distributions of kinematic variables [3(|. We describe the method in more detail in Sec. IIV Al 

The true 7 plus misidentified b background is estimated by first selecting events that pass all cuts except the 
requirement of a 6-taggcd jet, resulting in 18128 events (see Tab.|T|. For each selected event we then apply the product 
of two weights: (i) the true-photon weight determined using the CES/CPR method, representing the probability that 
a photon candidate is a photon, and (ii) the heavy-flavor mistag [3lj rate, which depends on jet Et, jet r\, the number 
of tracks in the jet, the number of primary interactions found in the event, and the z position of the primary interaction 
with the highest scalar sum pr of tracks. The mistag parameterization is the same as that used in the measurement 
of the ti cross-section [28|. Because the CES/CPR method and the mistag parametrization provide event-by-event 
weights, we are able to determine the shapes of kinematic distributions as well as the number of events for this 
background. 

We estimate the 76 and 7c backgrounds by generating MC events using madgraph [32| for leading-order matrix 
element processes involving photons, b or c quarks, and additional partons. The samples for 7 + b + jets and 7 + c + 
jets are generated with 1 to 3 jets. These samples are then processed with PYTHIA [33| to incorporate parton showering 
and hadronization. We ensure that we do not double-count events due to the overlap between jets arising from matrix 
element partons and jets arising from initial and final state radiation (341 ]. We obtain the overall normalizations of 
these backgrounds by fitting the secondary vertex mass distribution of the tagged jets, m(SV), to templates formed 
from the mass distributions of the expected SM components. The normalization scheme is described in Sec. IIV SI 

A summary of the background contributions is given in Tab. [TTJ 
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A. Photon Backgrounds: The CES/CPR Method 

For photon candidates with Et < 35 GeV, we use the shape of the shower profile measured with the CES system to 
discriminate between true single photon events and diphoton final states from decays of mesons. We construct a x 2 
discriminant by comparing the measured shower profile with that measured in electron test beam data [30] . A single 
photon has an average probability of w 78% to satisfy the x 2 cut, while the background has an average probability of 
~ 30% to satisfy the x 2 cut, since the shower profile of the two near-by photons from a meson decay is measurably 
wider on average. 

Above 35 GeV, however, the two photons from meson decay coalesce and the discrimination power of the shower 
profile measurement is significantly reduced. In this Et range, we use hit rates in the CPR system to discriminate 
between single photons and diphotons from meson decays. A single photon will convert and leave a hit in the preshower 
detector with a probability of w 65%. Backgrounds that decay into two photons have a hit probability of « 85% 
because the probability that neither photon converts is lower than the probability that a single photon does not 
convert. 

The difference of probabilities between signal (single photons) and background (photon pairs) forms the basis of a 
statistical method which assigns each event a weight for being a true photon (termed true-photon weight), W 7 , as 
described in Ref. (3(j. The weight is defined as: 

yy ^candidate ^bkg ^\ 

£sig — Cbkg 

where e s j g and ebkg are the respective probabilities for a true photon and misidentificd photon to satisfy a CES x 2 
cut or to leave a hit in the CPR, and (^candidate is either zero or one depending on whether the observed candidate 
satisfies these CES and CPR conditions. The values of e S i g and ebkg are determined using control data samples (3f| 
and are parametrized as a function of the energy of the photon candidate, the angle of incidence, and the number of 
primary interactions found in the event. The misidentified-photon weight is 1 — W 7 . We estimate the misidentificd 
7 background by summing up the misidentified-photon weights of the 617 candidate events (767^^) to obtain 115 
events with a statistical uncertainty of 49 events. We estimate the true 7, misidentificd b background by summing up 
the products of true-photon weights and heavy-flavor mistag rates of the 18128 events before the 6-tagging selection 
(77j'E T ) to obtain 141 events with a statistical uncertainty of 6 events. The calculation of the systematic uncertainty 
on these expectations is given in Sec. |Vj 

B. Heavy-Flavor Normalization 

The invariant mass of the tracks that form a secondary vertex can be used to discriminate between the bottom, 
charm, and light partons that compose a sample. We use this discriminating variable to normalize the contributions 
of the 76 and 7c backgrounds by fitting the secondary vertex mass distribution. 

The fitting technique utilizes templates of the distribution of the secondary vertex mass arising from the three 
primary sources expected to contribute to the observed distribution: bottom quarks, charm quarks, and light quarks 
or gluons. These templates are obtained from Monte Carlo samples containing final state photons [36|. The discrim- 
inati ng p ower of the secondary vertex mass is shown in Fig. [TJ in which the three templates are normalized to unit 
area [37J. The sum of the fractions of the three components is constrained to unity in the fit, /& + f c + /light = 1- 

This technique can be used to determine the number of events containing a real photon and real heavy flavor in any 
sample. We first subtract the contribution due to misidentificd photons by applying the CES/CPR method to obtain 
the number of misidentificd photon events. We then estimate the fraction of heavy flavor in events with a misidentificd 
photon by fitting the secondary vertex mass distribution in a sample enriched with jets faking photons, referred to 
as the sideband photon sample [38j . We then subtract the number of events containing a misidentified photon and 
heavy flavor from the number of events obtained from the standard photon sample fit to obtain the number of 76 and 
jc events. 

In principle this technique could be directly applied to the search region to obtain the number of 76 events. However, 
this would not give us the SM expectation, as the contribution from any new process making such events would be 
counted [3!| . Instead, the expected 76 contribution is normalized by applying this technique to a control region with 
a much larger SM cross-section than that of the search region, and then extrapolating to the search region by using 
efficiencies derived from the 7 + 6 Monte Carlo. The final estimate for the number of 76 events in the search region 
is (search) = iV 7 f, (control) • e (control — > search). 

We define the control region as the 7 + 6-tag sample, where the only selection requirements are that there be at 
least one photon with \t}\ < 1.1 and E T > 25 GeV and one SECVTX-tagged jet having |r?| < 2 and E T > 15 GeV. The 
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FIG. 1: Templates of the invariant mass of all tracks in a secondary vertex arising from bottom quarks, charm quarks, and 
light quarks and gluons from Monte Carlo simulations, normalized to unit area. 



number of events in the search region is less than 1% of that in the control region, which contains 93 894 events. We 
obtain an efficiency of ^(control — * search) = 0.0123 ± 0.0025, defined as the fraction of 7 + b Monte Carlo events in 
the control region that survive the additional cuts of the search region. The uncertainty on the efficiency is due to 
the differences in jet multiplicities and E T distributions between data and the background prediction in the control 
region. 

Figure [2] shows the results of a maximal likelihood fit performed on the search and control region using the templates 
above to extract the fraction of &-jet and c-jet events. We estimate the number of 76 events by subtracting the 
misidentified photon plus b contribution from the control region and then multiplying by e (control — > search) to 
obtain 341 events with a statistical uncertainty of 18 events. The calculation of the systematic uncertainty on the 
number of events is given in Sec. |V] 

The 7c background is normalized by directly fitting the secondary vertex mass in the search region. We do not 
extrapolate the charm normalization from the control sample because the uncertainties on the matching scheme for 
charm quarks are large (40| and therefore the extrapolation efficiency would have large uncertainties. After subtracting 
the misidentified photon plus charm contribution we obtain an estimate of 9 7c events with a statistical uncertainty 
of 52 events. 

Note that because the charm background is measured in the search region, this search is not sensitive to anomalous 
charm production. It is, however, sensitive to anomalous production of the 767 E T final state because we use 76 Monte 
Carlo processes to obtain the efficiency e(control — > search). 



V. SYSTEMATIC UNCERTAINTIES 



Sources of systematic uncertainty on the number of predicted events arise from: a) the uncertainty on the true- 
photon weights W 7 in the CES/CPR method, b) the uncertainty on the heavy-flavor mistag prediction, and c) the 
uncertainty on the template shapes used in the secondary vertex mass fit. 

The systematic uncertainty on W 7 in the CES/CPR method arises from uncertainties on the CES x 2 efficiencies 
and the CPR hit rates for photons and backgrounds [e s i g and 6bkg in Eq. {!])]. The largest uncertainty on the CES 
X 2 efficiencies is due to the gain saturation in the CES detector [3fJ (10% on W 7 ). The largest uncertainty on the 
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FIG. 2: The secondary vertex mass fit in events containing standard photons, for the control sample (left) and the search 
sample (right). Note the uncertainties on f c and ft, are correlated and purely statistical. 



CPR hit rates is due to the modeling of the hit rates of tt° and r/° (5-35% on W 7 ). On average, W 7 has a relative 
systematic uncertainty of 11%. 

The uncertainty on the heavy-flavor mistag prediction comes from: the finite size of data samples used for parame- 
terization of the mistag rates (10%), variations between different data-taking periods (6.5%), and the uncertainty on 
a scale factor which takes into account the contribution of misidentified 6-tags from long-lived hadrons (A , K®) and 
secondary particles due to interactions with detector material (10-15%). More detail may be found in Ref. |28f . 

We estimate the systematic uncertainty arising from the secondary vertex mass fitting procedure by varying the 
shapes of the templates that are used in the binned likelihood fit. The systematic effect of mismodeled tracking 
inefficiency in the Monte Carlo simulation is estimated by lowering the secondary vertex mass template mass scale 
by 3% (41J. We also refit the secondary vertex mass distributions with templates derived from Monte Carlo samples 
that have the E T > 25 GeV cut imposed on them, as this may change the relative fraction of semileptonic decays 
in the template samples and thereby alter the secondary vertex mass distribution. Because both of these sources 
of uncertainty affect the shape of the templates, we take the maximum variation observed as the systematic shift 
in normalization. We obtain a 12% uncertainty on the 6-fraction and a 48% uncertainty on the c-fraction from this 
estimate. 

The numerical values of the systematic uncertainties are presented in Tab. [II] in Sec. I VII below. The CES/CPR 
method contributes 13% of the systematic uncertainty on the total amount of background while the mistag 
parametrization and secondary vertex mass fit contribute 24% and 63% respectively. The calculation of the total 
systematic uncertainty takes into account correlations among the different sources of backgrounds to the jbj^r sig- 
nature. Because the CES/CPR method is used to estimate the contribution of all four background categories defined 
in Sec. lIVi we apply the CES/CPR systematic variations to all backgrounds simultaneously when calculating the final 
CES/CPR uncertainty on the total background prediction. All other sources of systematic uncertainty are combined 
as uncorrclated uncertainties. 



VI. RESULTS 



We proceed to test the SM in the jbj^T signature in three ways: comparing predicted event counts, looking 
for anomalous kinematic behavior, and counting additional objects in the events, as might be expected from the 
production of new heavy states with extended decay chains. We also go beyond the Run I measurement criteria by 
increasing the requirement on missing transverse energy to 50 GeV, reducing the expected background contribution 
by a factor of ss 20, and thereby enhancing the sensitivity to new processes. The three tests are described in the 
sections below. 
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A. Comparing Predicted Event Counts 

Tabic |TT] summarizes the background sources with associated statistical and systematic uncertainties. The total 
background prediction is 

N(BG) = 607 ± 74 (stat.) ± 86 (syst.), (2) 

where the first uncertainty is statistical and the second systematic. The observed number of events is 617, consistent 
with the background predictions. 

TABLE II: The numbers of predicted events from background sources. The two uncertainties in each row are statistical and 
systematic, respectively. Note that the total systematic uncertainty is less than the largest individual contribution due to an 
anti-correlation of the CES/CPR uncertainties between the components. 



Background Source 


Expected Events 




Stat 




Sys 


Misidentified 7 


115 


± 


49 


± 


54 


True 7, misidentified b 


141 


± 


6 


± 


30 


") b 


341 


± 


18 


± 


91 


7c 


9 


± 


52 


± 


14 


Total 


607 


± 


74 


± 


86 



B. Object Kinematics 

Wc examine three different types of distributions for anomalous shape discrepancies with respect to the background 
prediction: the kinematics of individual objects in the event such as jets and photons, global features of the event 
such as J5 T , and the invariant masses of the combinations of objects. 

The distributions of the transverse energy of the photon, 6-jet, and 2 nd jet are shown in Figs. [3l [4J and [5] respectively. 
The distributions of J5 T , TV(jets), and Ht, where Ht is the scalar sum of the transverse momenta of the photon, 
all jets in the event, and J5 T are shown in Figs. O and [5] respectively. The J£ T distribution is shown before the 
application of the J5 T > 25 GeV cut but after the application of all other selections. The distributions of M(-fb), 
M(bj), M(~fbj), M T (7p T ), and M T (bjfi T ) arc shown in Figs. M EH CLU HH and [T3l respectively. The transverse mass 
Mt is calculated with the transverse components of object momenta: 

M T = ^(I]^) 2 -^) 2 - (S^) 2 , (3) 

where E T , p l x , and p l y are the transverse energy and x and y components of the momentum of object i (which could 
be a photon, 6-quark jet, jet, or missing energy). Note that the binning for all distributions is such that there are no 
overflows. 

We test the consistency between the observed shapes of kinematic distributions and the shape predicted by the 
background expectation by running pseudocxperiments for each distribution studied and calculating the Kolmogorov- 
Smirnov (KS) distance for each pscudocxpcrimcnt. The probability that a random sampling of the estimated back- 
ground distribution would give a higher KS distance than the observed data distribution, referred to as the "KS 
p- value" , is obtained for each kinematic variable studied by integrating the tail of the distribution of KS distances. 
We obtain a range of KS p-values between 7.0% and 99.8%, indicating that the kinematic distributions observed are 
consistent with background expectations. 

C. Effect of Additional Selections 

Wc further investigate the existence of possible anomalies in the 767 ,E T X final state by making additional selections 
and comparing the number of observed events to the background predictions. We chose criteria based on expected SM 
distributions and selections used previously in the search of Ref [T(| • The additional selections we make are J5 T > 50 
GeV, TV (jets) > 3, p T (-y) > 50 GeV, H T > 200 GeV, E T (b) > 50 GeV, and A0(jet, £ T ) > 0.5. 

Table UTTI summarizes the effects of the additional selections. We apply these in two different ways: one at a time 
independently of all other additional selections, and after the application of the J£ T > 50 GeV selection. No anomalous 
excess of events is observed. 




E T ( Y ) [GeV] 

FIG. 3: The distribution in photon Et observed (points) and from backgrounds (histogram). The KS p- value is 63.7%. 




E T (b) [GeV] 

FIG. 4: The distribution in 6-jet Et observed (points) and from backgrounds (histograms). The KS p- value is 59.7%. 



Finally, exotic particles with cascade decays, X — > jY — > bj, may form a cluster in the scatter plot of M(^bj) vs. 
M{bj). In Fig. [HI we compare the observed distribution to that from the estimated background; we do not see any 
evidence of an anomaly. 
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E T (j) [GeV] 

FIG. 5: The distribution in untagged jet Et observed (points) and from backgrounds (histogram). The KS p-value is 10.4%. 




FIG. 6: The distribution in missing transverse energy observed (points) and from backgrounds (histogram). The KS p-value is 
7.0%. 



VII. CONCLUSIONS 



We have searched for the anomalous production of events containing a photon, two jets, one which is identified as 
originating from a b quark, and missing transverse energy. The number of events observed in data is consistent with 
the number of expected background events. No significant excess of events with respect to the background prediction 
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FIG. 7: The distribution in jet multiplicity observed (points) and from backgrounds (histogram) in logarithmic scale. The KS 
p-value is 19.0%. 



TABLE III: The number of events observed and the predicted background for additional independent selections. The first 
uncertainty in the observed columns is statistical and the second is systematic. 



Selection No additional cuts With E T > 50 GeV 

Observed Predicted Observed Predicted 

E T > 50 GeV 28 30 ± 10 ± 5 

iV(jets) > 3 321 329 ± 46 ± 46 15 17 ± 7 ± 3 

p T (j) > 50 GeV 257 247 ± 42 ± 39 16 21 ± 8 ± 5 

H T > 200 GeV 304 322 ± 45 ± 46 25 28 ± 9 ± 5 

E T {b) > 50 GeV 286 310 ± 43 ± 44 18 22 ± 8 ± 6 

A^(jet, E T ) > 0.5 343 368 ± 47 ± 49 15 16 ± 8 ± 4 



is observed in any of the kinematic distributions studied. The shapes of these distributions are consistent with SM 
expectations. Furthermore, we do not see any anomalous production of events after applying additional selections. 
We conclude that the 2.0 ±0.1 fb -1 7 + b + j + E T + X sample is consistent with SM background expectations. 
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E T (y) + Z E T (j) + £ T [GeV] 

FIG. 8: The distribution of the scalar sum of the transverse momenta of the 7, all jets in the event, and Fj t observed (points) 
and from backgrounds (histogram). The KS p- value is 99.7%. 



and Programa Consolider-Ingenio 2010, Spain; the Slovak R&D Agency; and the Academy of Finland. 



15 



. Data, J/, dt = 2.0 fb" 1 

I I Fake y, Real+Fake b 
I I Real y, Fake b 

□ yb 

I" 7777 ! Background Uncertainty 




M(yb) [GeV/c ] 



FIG. 9: The distribution of the mass of the photon + b jet observed (points) and from backgrounds (histogram). The KS 
p- value is 62.0%. 



. Data, J/, dt = 2.0 fb" 1 

I I Fake y, Real+Fake b 
I I Real y, Fake b 

□ yb 

I I Background Uncertainty 




100 200 



700 

M(bj) [GeV/c 2 ] 



FIG. 10: The distribution of the dijet mass observed (points) and from backgrounds (histogram ). The KS p-value is 9.8%. 
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FIG. 11: The distribution of the invariant mass of the 7, fo-jet, and 2 nd jet observed (points) and from backgrounds (histogram) 
The KS p-value is 99.8%. 
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FIG. 12: The distribution of the transverse mass of the photon + JC T observed (points) and from backgrounds (histogram). 
The KS p-value is 96.8%. 
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FIG. 13: The distribution of the transverse mass of the &-jet, 2 nd jet, and missing transverse energy observed (points) and 
from backgrounds (histogram). The KS p-value is 21.1%. 
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FIG. 14: M(bj) versus M("/bj) for the events which satisfy the selections in Tab. [I] observed (big dots) and expected (small 
dots). 
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